A quantitative model for IqR product in d-wave Josephson junctions 
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We study theoretically the Josephson effect in d-wave superconductor / diffusive normal metal 
/insulator/ diffusive normal metal/ d-wave superconductor (D/DN/I/DN/D) junctions. This model 
is aimed to describe practical junctions in high- Tc cuprate superconductors, in which the product of 
the critical Josephson current (7c) and the normal state resistance (R) (the so-called IcR product) 
is very small compared to the prediction of the standard theory. We show that the IcR product 
in D/DN/I/DN/D junctions can be much smaller than that in d-wave superconductor / insulator 
/ d-wave superconductor junctions and formulate the conditions necessary to achieve large IcR 
product in D/DN/I/DN/D junctions. The proposed theory describes the behavior of IcR products 
quantitatively in high- To cuprate junctions. 



Josephson effect in high Tq superconductors has at- 
tracted much attention-^ because of its potential ap- 
plications in future technologies^. In particular, ap- 
plications in electronics, such as the single-flux quan- 
tum devices are extremely promising, since the oper- 
ating frequency is proportional to the product of the 
critical Josephson current (Jc) and the normal state re- 
sistance (R) (the so-called IcR product) which is ap- 
proximately proportional to the superconductivity criti- 
cal temperature 4 . However, almost all experimental data 
for high-Tc Josephson junctions fabricated so far have 
shown IqR values much smaller than those predicted by 
the standard theory, irrespective of what kind of junc- 
tions they are££. This strongly suggests that the in- 
terfaces in the cuprates are intrinsically pair breaking. 
Several models were proposed which treat this issue^&iL. 
However, the situation was not described so far in a quan- 
titative manner. 

One possibility addressed in the present paper and not 
investigated before is that superconductivity is destroyed 
near the interface in d-wave superconductor / insulator / 
d-wave superconductor (DID) junctions, where diffusive 
normal metal (DN) regions are induced. Thus, DID junc- 
tions turn into d-wave superconductor / diffusive normal 
metal /insulator/ diffusive normal metal/ d-wave super- 
conductor (D/DN/I/DN/D) junctions (see Fig. In 
these junctions, IcR product can be much smaller than 
that in DID junctions. In this paper we will explore 
this possibility and provide a quantitative model which 
is compared to the experimental data. 

The Josephson effect is a phase-sensitive phenomenon 
and thus depends strongly on a superconducting pairing 
symmetryii2ii£ . In DID junctions, nonmonotonic tem- 
perature dependence of critical curren t 11 i 12 ' 13 i 14 ' 15 occurs 
due to the formation of midgap Andreev resonant states 
(MARS) at the interface^. The MARS stem from sign 
change of pair potentials of d-wave superconductors 1 ^. 
It was also predicted that MARS strongly enhance the 



Josephson current at low temperatures^- 2 -. On the other 
hand, in Josephson junctions with DN, the role of the 
MARS change. 

In superconductor / diffusive normal metal / supercon- 
ductor (S/DN/S) junctions Cooper pairs penetrate into 
the DN as a result of the proximity effect, providing the 
Josephson couplin g 18 i 19 ' 20 i 21 ' 22 . Scattering of electrons 
by impurities in the DN layer makes superconducting co- 
herence length shorter and thus suppresses the Joseph- 
son current. In D/DN/D junctions, the Josephson cur- 
rent is suppressed by the MARS 24 ' 25 ' 26 ' 27 , in contrast to 
DID junctions, because MARS compete with proximity 
effect^^. Therefore I C R product in D/DN/I/DN/D 
junctions can be much smaller than that in DID junc- 
tions. 

In the present paper, we calculate Josephson current 
in D/DN/I/DN/D junctions as a model of the actual 
DID (e.g., grain boundary) junctions. We show that 
IcR product in D/DN/I/DN/D junctions can be much 
smaller than that in DID junctions, and clarify the con- 
ditions with which the IcR product is most enhanced in 
D/DN/I/DN/D junctions. Our theory can explain the 
above mentioned general trend of the high-Tc Josephson 
junctions quantitatively, in contrast to previous theoret- 
ical models of high-Tc cuprate junctions. The obtained 
results may provide useful information for fabrication of 
high-Tc Josephson junctions. 

Let us formulate the model for a D/DN/I/DN/D junc- 
tion. We assume that the DN layer has a length L 
much larger than the mean free path and is character- 
ized by the resistance Rd- The DN/D interfaces lo- 
cated at x — ±L have the resistance R' b , while the 
DN/I interface at x = has the resistance We 
model infinitely narrow insulating barriers by the delta 
function U(x) = H'5{x + L) + H8{x) + H'5{x - L). 
The resulting transparencies of the interfaces T m and 
T' m are given by T m = 4 cos 2 </>/(4cos 2 + Z 2 ) and 
T' m = 4cos 2 0/(4 cos 2 <j> + Z' 2 ), where Z = 2H/v F and 
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Z' = 2H'/vp are dimensionless constants and vp is 
Fermi velocity, where <f> is the injection angle measured 
from the interface normal. In the following we assume 
Z 3> 1. The schematic illustration of the model is 
shown in Fig. [TJ The pair potential along the quasi- 
particle trajectory with the injection angle (f> is given by 
A L = Acos[2(^-a)]exp(-z*) and A R = A cos[2(0-/3)] 
for the left and the right superconductors, respectively. 
Here *S? is the phase difference across the junction, a and 
(3 denote the angles between the normal to the interface 
and the crystal axes of the left and right d-wave super- 
conductors, respectively. The lobe direction of the pair 
potential and the direction of the crystal axis are chosen 
to be the same. 
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FIG. 1: (color online) Schematic illustration of the model for 
the D/DN/I/DN/D junction. 



We parameterize the quasiclassical Green's functions 
G and F with a function $.i2a±2: 



G^ 



y/u 2 + p!. 



- F 



y/tU 2 + $ w $*_ u 



(1) 



where u> is the Matsubara frequency. In the DN layers 
the Green's functions satisfy the Usadel equation 23 



where T is temperature and R = 2R d + Rb + 2R' b is the 
normal state resistance of the junction. In the following 
we focus on the IcR value as a function of temperature 
and clarify the cases when IcR is enhanced. Below A(0) 
denotes the value of A at zero temperature. Note that it 
is realistic to choose small magnitude of Z' and R' b , and 
large Thouless energy because thin DN regions could be 
naturally formed due to the degradation of superconduc- 
tivity near the interface. 
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where £ = y/ D /2irTc is the coherence length, D is the 
diffusion constant and Tq is the transition temperature 
of superconducting electrodes. To solve the Usadel equa- 
tion, we apply the generalized boundary conditions de- 
rived in Ref j 24 i 25 at x = ±L and the boundary conditions 
in Ref. 21 at x = 0. 

The Josephson current is given by 
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FIG. 2: (color online) I C R value for Z' = 1, R d /Rb = 0.1 and 
(a,/?) = (0,0). 

In FigEJwe show I C R value for Z' = 1, R d /R b = 0.1 
and (a,/3) = (0,0) with various E Th /A(0) and R d /R' b . 
IcR increases with Etii/A(0) and Rd/R' b because prox- 
imity effect is enhanced. As Eth increases, the magni- 
tude of the gradient becomes small. 

Figure [3] shows I C R value for E Th /A(0) = 1, R d /R b = 
0.1 and (a, j3) = (0, 0) with various Z' and R d /R' b . As Z' 
increases, the magnitude of the gradient becomes large. 
The peculiar effect is that IcR increases with Z', indi- 
cating that proximity effect is enhanced by the increase 



3 




FIG. 3: (color online) Ic R value for E Th /A(0) = 1, R d /R b = 
0.1 and (a, 0) = (0,0). 

of Z' . This stems from the sign change of the pair 
potentia l 24 i 25 . For the case of d-wave symmetry with 
a = = 0, injection angles of a quasiparticle can be 
separated into two regions: cf) + = {0|O < \<f>\ < ir/i} and 
= {(/)\tt/4 < \<j>\ < 7r/2}. The signs of pair potential 
for (f) + and that for cj>_ are opposite. As a result, the 
sign change of pair potentials suppresses the proximity 
effect in the DN and hence Josephson currents. As Z' 
increases, the contribution from <f> + dominates over that 
from 0_. Therefore IcR increases with Z' . 

In Fig. Iwe plot I C R value for E Th /A(0) = 1, Z' = 
1 and (a,0) = (0,0) with various Rd/Rb and Rd/R' b - 
IcR increases with Rd/Rb due to the enhancement of 
the proximity effect. 

Figure [5] displays I C R value for E Th /A(0) = 0.1, 
Z' = 0.1, Rd/Rb = 0.1 and R d /R' b = 10 with various 
a and 0. The formation of MARS suppresses the prox- 
imity effect. Therefore IcR decreases with the increase 
of a and / 3 24 i 25 i 26 i 27 , In the actual junctions, there is in- 
evitable roughness at the interface and hence the effective 
values of a and at the interface become random even if 
junctions with a — = are fabricated. This provides 
the mechanism of suppression of the IcR product. 

Finally we compare the present theory with the ex- 
perimental data from Ref£ and with the theory for DID 
junctions by Tanaka and Kashiwaya (TK)^ 2 -. The tem- 
perature dependencies of Ic are plotted in Fig. [6] taking 
a = = and R = 0.375S1 for theoretical plots. We 
choose E Th /A(0) = 3, Z' = 0.1, R d /Rb = 0.01, and 
Rd/R'b = 100 in the present theory, and the barrier pa- 
rameter Z = 10 in the TK theory. As shown in this 
figure, the present theory can explain the experimental 
results quantitatively, while the discrepancy between the 
TK theory and the data is rather strong, about an order 
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FIG. 4: (color online) I C R value for E Th /A(0) = 1, Z' = 1 
and (a, 13) = (0,0). 
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FIG. 5: (color online) I C R value for E Th /A(0) = 0.1, Z' = 
0.1, R d /R b = 0.1 and R d /R' b = 10. 



of magnitude. Note that in the TK theory the Ic is not 
sensitive to the choice of Z parameter. To estimate the 
realistic size of the DN region, we can take A(0) = lOmeV 
and D = 10 _3 m 2 /s, and then obtain the length of the 
DN region L = 4.7nm. 

In summary, we have studied the Josephson current in 
D/DN/I/DN/D junctions as a model of high Tc super- 
conductor junctions. We have shown that the IcR prod- 
uct in D/DN/I/DN/D junctions can be much smaller 
than that in DID junctions and have found the conditions 
when the IcR in D/DN/I/DN/D junctions is largest. 
The requirements for the large magnitude of IcR prod- 
uct are: no roughness at the interfaces, large magnitudes 
ofZ', R d /Rb, R d /R' b and E Th , and (a,0) = (0,0). Note 
that small magnitude of Z' and R' b , and large Exh are 
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FIG. 6: (color online) Comparison between the present theory 
(solid line), experimental data 6 (dotted line) and TK theory 12 
(broken line). Ic is plotted as a function of temperature, 
taking R — 0.375f2 and a = (3 = for theoretical plots. 
We choose _E Th /A(0) = 3, Z' = 0.1, R d /R b = 0.01, and 
Rd/R'b = 100 in the present theory, and Z — 10 in TK theory. 



realistic for naturally formed DN layers, hence the only 
tunable parameter is Rh. Our theory can explain the ex- 
perimental results on the quantitative level, in contrast 
to the previous idealized treatment of DID junctions. 
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